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ABSTRACT 


[G],[Gs] matrices defined In Eq. (6) 


The effects of structural coupling on mlstuned 
cascade flutter and response are analytically Investi- 
gated using an extended typical section model. Previ- 
ous work using two degree of freedom per blade typical 
section monels has included only aerodynamic coupling. 
The present work extends this model to include both 
structural and aerodynamic coupling between the blades. 
The model assumes that the structurally coupled system 
natural modes have been determined and can be repre- 
sented in the form of N bending and N torsional un- 
coupled modes for each blade, where N is the number of 
blades and, hence, is only valid for blade dominated 
motion. The aerodynamic loads are calculated by using 
two-dimensional unsteady cascade theories in the sub- 
sonic and supersonic flow regimes. The results show 
that the addition of structural coupling can affect 
both the aeroelastic stability and frequency. The 
stability is significantly affected only when the 
system is mistuned. The resonant frequencies can be 
significantly changed by structural coupling in both 
tuned and mlstuned systems, however, the peak response 
is significantly affected only in the latter. 
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K hsr* K asr 



defined In Eq. (6) 

bending deflection of blade In r** 1 
mode of tuned cascade 

bending deflection of s** 1 blade 

mass moment of inertia of s th blade 
about elastic axis per unit span 

square root of minus one, also a 
summation Index 

bending and torsional stiffness of 
s th blade In r™ mode, respectively 

reduced frequency based on semi chord 

lift due to motion of s th blade 
per unit span, positive up 

lift due to wakes of s th blade per 
unit span, positive up 


NOMENCLATURE 

[A] aerodynamic matrix due to motion 

(AD) aerodynamic matrix due to wake In- 

duced flow 

a elastic axis position 

b semichord 

c chord 


^hhr»^har 


^hr* 'a«r 
^hr* ^war 


nondimens ional lift coefficients due 
to bending and torsional motions, 
respectively, in the r th mode 

nondlmensional moment coefficients 
due to bending and torsional motion, 
respectively, ir. the r^ h mode 

nondlmensional lift and moment coef- 
ficients, respectively, due to wakes 
in the r* h mode 


ro],[D s ] matrices defined in Eq. (3) 

[E] transformation matrix 

e base for natural logarithm 
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Mach number 

moment about the.elastic axis due to 
motion of the s tn blade per unit 
span, positive nose up 

moment of s th blade per unit span 
about the elastic axis due to wake, 
positive nose up 
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mass per unit span of the s th blade 

number of blades 

matrix defined in Eq. (10) 

matrix defined in Eq. (7) 

Integer specifying the nodal diame- 
ter or Interblade phase angle mode 
of the tuned rotor 

radius of gyration of s th blade, 
nondimensionallzed with respect to b 

integer specifying blade number, 
also blade . 

static mass moment of the s th blade 
per unit span about elastic axis, 
positive when center of gravity is 
aft of the elastic a^is 

time 

freestream velocity relative to the 
blade 

physical blade displacements 
rectangular coodlnate axes 

static unbalance of 


di ..acements in nodal diameter modes 

amplitude of torsional motion of 
s th blade, positive clockwise 

torsional amplitude of each blade of 
the tuned rotor 

amplitude of torsional deflection of 
a blade in the r t ' 1 mode of a tuned 
cascade 

interblade phase angle, 2»r/N 

nondimensional eigenvalue, (« 0 /“) 2 

defined in Eq. (14) 

defined in Eq. (4) 

damping ratios of the s^ blade in 
bending and torsion in r th mode, 
respectively 

mass ratio of the s th blaoe 

real part of eigenvalue, defined in 
Eq. (12) 

imaginary part of eigenvalue, de- 
fined in Eq. (12) 

stagger angle 

frequency 

reference frequency 


nondimensional 
s l Made 


“hs 


“os 

V’OTs 

“hsr 

bending natural frequency, dependent 
on r 

“osr 

torsional natural frequency, depen- 
dent on r 

“hr 

tuned bending frequency, dependent 
on r 

“or 

tuned torsional frequency, dependent 
on r 

[ H 

Inverse of a matrtx 


I. INTRODUCTION 

In Refs. (1) and (2) the authors have Investigated 
the effects of mistuning on cascade flutter and re- 
sponse by assuming that the Individual blades in the 
cascade are coupled only aerodynamical 1y. A typical 
section model, in which each blade has two degrees of 
freedom, one bending *nd one torsional, was used. In 
all bladed-disc assemblies the blades are structurally 
coupled through the disc and often through a variety 
of other connecting parts such as shrouds, dampers, 
and lacing wires. In the case where there Is only weak 
structural coupling between the blades (e.g., a rela- 
tively stiff disc and flexible blades), the effect of 
the structural coupling on individual blade frequencies 
is negligible and the method of analysis used in the 
previous publications Is adequate. To handle cases 
where there is significant structural coupling, the 
typical section model must be refined. The purposes 
of this paper are to present such a refined model to 
account for structural coupling and to study Its ef- 
fects on cascade flutter and aeroelastlc response of a 
mistuned bladed-disc. The work presented is limited 
to bladed-discs In which the motions are dominated by 
the blades and is based on a portion of the research 
described in Ref. (3). 

The only known"published work on mistuning in- 
cluding both aerodynamic and structural coupling is 
given In Ref. (4), In which each blade is allowed a 
single torsionaT degree of freedom and a subsonic aero- 
dynamic theory is used to investigate flutter. In the 
present paper the model is capable of considering cou- 
pled bending and torsional motion and both subsonic 
(Ref. (5)) and supersonic (Ref. (6)) two dimensional 
unsteady cascade aerodynamic theories are used to in- 
vestigate flutter and aeroelastlc response. Since the 
effects of bending-torsion coupling were examined in 
Refs. (1) and (2), the parametric studies are limited 
to predominantly torsional motion with structural cou- 
pling. In addition to the references mentioned above, 
some other publications dealing with mistuning are (1) 
structural coupling without aerodynamics (Refs. (7) to 
(10)) and (2) flutter and response using single degree 
o7~freedom blade models with aerodynamic coupling only 
(Refs. (_1_1) (14)). Tuned systems with two degrees 

of freedom per bTade have been considered in Refs. (15) 
to (20). In Ref. (21) the authors have considered a 
more complex beam model which is capable of repre- 
senting a twisted nonuniform blade. 

II. THEORY 


o 


2 


It is well known that the natural frequencies of 
annular plates depend upon the number of nodal diame- 
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ters. In many publications (e.g., Refs. (2) to (10)) 
this behavior has also been shown to be true for 
bladed-disc assemblies. For tuned systems each single- 
blade mode, such as the first bending mode, becomes a 
family of inodes, one for each nodal diameter pattern. 

In general, the frequencies for the first bending fam- 
ily of modes can vary considerably and depend monoton- 
Ically on the number of nodal diameters. However, for 
torsional motion the frequencies are relatively Insen- 
sitive to nodal diameter and do not increase monoton- 
Ically. For a more detailed discussion of this behav- 
ior see Ref. (3). 

In general, a complete bladed-disc model must in- 
clude blade models, a disc model, and the appropriate 
continuity conditions. Simpler approaches can be used 
for the two limiting conditions of blade and disc dom- 
inated motipn. For disc dominated motion, the blades 
can be considered to be lumped masses with no effect 
on the system stiffness. For blade dominated motion, 
the disc can be considered to be a spring support with 
no effect on the system mass. This latter approach is 
used herein and, hence, is only valid for blade domi- 
nated motion. The authors are developing a complete 
bladed-disc model which will be described in a future 
publ Ication. 

The approach used In formulating the mathematical 
model Is to assume that the structurally coupled sys- 
tem modes (In a vacuum) have been determined and can 
be represented In the form of N bending and N tor- 
sional uncoupled modes for each blade. That is, each 
blade will appear to have 2N natural frequencies: 
one bending and one torsion for each nodal diameter 
pattern. This description requires that the disc (and/ 
or shroud ring) De a perfect .body of revolution with 
all of the mistuning effects lumped in the Individual 
blades. Since the number of nodal diameters, r, is 
directly related to the interblade phase angle by 

8 r - 2»r/N (1) 

each interblade phase angle mode is associated with a 
given nodal diameter. The parameter, 8 r , is used in 
the "traveling wave approach" of describing blade mo- 
tion. In this approach the motion of the airfoils in 
each mode of the tuned cascade is assumed to be simple 
harmonic with a constant phase angle between adjacent 
blades. This approach is entirely equivalent to the 
"standing wave approach" in which the motion Is char- 
acterized by modes with nodal diameters fixed in the 
disc. For further discussion of the standing and 
traveling wave approaches see Ref. (20). To include 
structural coupling, the equations oT” motion given In 
Ref. (c) have to be modified to reflect the fact that 
the blade frequencies in a vacuum depend on the number 
of nodal dia..«ters. This is accomplished by assuming 
that only the stiffness terms are functions of the 
number of nodal diameters. Hence, the natural fre- 
quencies in a vacuum are given an additional subscript 
to reflect this dependency. Also, the force (or mo- 
ment) associated with the stiffness In the equations 
of motion is assumed to be proportional to the dis- 
placement in the r th nodal diameter mode, and not to 
the physical blade deflection. As a result, the 
stiffness terms must be multiplied by the modal dis- 
placement and not the blade physical displacement. 

The term "modal displacement" is used since the physi- 
cal displacements are expanded in terms of the nodal 
diameter modes. Also, the structural damping associa- 
ted with the bending and co-sional motions is also 
allowed to be a function of the number of nodal 
di ameters. 

as in Refs. (_1) and (2), each blade is modeled as 
a two degree of freedom oscillator in which the plung- 
ing and pitching motions are i nert i a 1 1 y coupled (see 


Fig. 1). Including the nodal diameter dependence de- 
scribed above, the equations of motion of the s th 
blade are 



These are Identical to Eq. (4) of Ref. (2) with the 
following exceptions: the stiffness termc are now ex- 
panded to Include the stiffness associated with each 
nodal diameter mode and the physical displacement terms 
are replaced with the modal displacements. Note that 
the real parts of the modal stiffnesses, K^sr and 
K oSr , are expressed In the form of the mass times 
square of the frequency. The aerodynamic forces are 
the same as those given In Ref. (2) and are also de- 
pendent on the nodal diameter. Tlondlmenslonallzlng 
Eq. (2) results In 



N-l 


Z 


r«0 


(1 * 


^ + 21; hsr^ u s r as Y asr' r 


rh ar /b> » 


2»1rs 


> e 


N-l 


z 


r-0 


'hhr 

*har 

Jahr 

1 r 

aar J 


<V 6 ' 



-( 3 ) 


where 


3 



ORIGINAL PAGE IS 
OF POOR QUALITY 


Y hsr " “hsr /u o 


' asr 


"asr'-o 


(4) 


Y « (a» 0 /«i»)^ 

The other nondlmenslonal variables are defined In the 
nomenclature and described In Ref. (1). Writing these 
equations for all blades and arranging them In matrix 
form results In 

[D]{X> ♦ y [Q]<Y> - [G][E][A]{Y> ♦ [GJ[E]{AD.> (5) 

where 
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[0,3 
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The complex stiffness and mass matrices are repre- 
sented by [G] and [D], respectively. The (X) and 
<Y } vectors represent the physical blade displace- 
ments and the modal displacements, respectively. The 
[El, [A], and {AD} matrices are defined In Ref. 

(2) and are not repeated here for brevity. Using the 
relationship 

(X) - [E]{ V} ( 8 ) 

to remove the physical blade displacements from the 
system of equations results In 

[P]<Y> - y[Q]{Y) - CG][E]{AD} (9) 

where 

[P] - [G][E][A] ♦ [DJ[E] (10) 

The aeroelastlc stability of the system, can be deter- 
mined from the eigenvalues, y's, of the generalized 
problem 

[fP] - y[Q]](Y} - 0 vll) 

The relation between the eigenvalue, y. and the com- 
plex frequency is 


i- . L_ „ 7 * 17 ( 12 ) 

“o V7 


Therefore, y is a measure of the damping or stability 
of the mode and v is a measure of the frequency. For 

the discussion of results to follow, a mode is consi- 
dered to be stabilized when M becomes more negative 
and destabilized when It becomes more positive. The 
aeroelastlc response can be derived from Eq. (9) and Is 

<Y> - [y[Q] - [P]]" l [G][E)[AD] (13) 

II. RESULTS ANO DISCUSSION 

A computer program was written to calculate the 
stability and response of a mlstuned bladed-disc In- 
cluding both aerodynamic and structural coupling of 
the blades using the formulation presented in the last 
section. In this program it is possible to consider 
arbitrary mistuning of the blade uncoupled bending and 
torsional frequencies and damping coefficients. The 
quantities Yhsr an<} Y a sr» see Eq . ( 4) , must be 
supplied as Input to this formulation. The approach 
used herein is to calculate the frequencies of the 
tuned rotor for each nodal diameter up to N/2 (N/2-1 
for an odd number of blades). 


4 
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Y hr " “hr /u o * Y hN-r 
Y ar " “ar /u o " Y aN-r 


(14) 


For example, for a 24 bladed rotor In vacuum the fre- 
quency of the r - 5 mode Is the same as that of the 
r - 19 mode since the structural properties cannot 
distinguish between a forward and backward traveling 
wave of the same noual diameter. There are a number 
of methods available to calculate the quantities 
vhsr and YasrS the most general being the finite 
element method. However, It Is sufficient for the 
purposes of this paper to assume that the mlstuned 
system frequencies are given by 


Y hsr " T hr(“hs / “o ) 
Y asr * Y ar ( “aS / “o ) 


(15) 


in which and <ii aS are the blade alone fre- 

quencies. This has been shown In Ref. (22) to be a 
valid assumption for small deviations in u^s and 
u> aS when the bladed-dlsc motion Is blade dominated. 

In the following sections two different rotors are 
analyzed: (1) a twelve bladed rotor representative of 
a set of blades connected by midspan shrouds (consi- 
dered In Ref. (4)) and (2) a more realistic unshrouded 
bladed disk representative of an advanced fan stage. 
Both alternating and random torsional frequency mis- 
tuning are considered. Since the effects of bending- 
torsion coupling were examined In Refs. (1) and ( 2 ), 
the emphasis in this paper is on the torsional mo- 
tion. As a result the predominantly bending modes are 
not shown in the results to follow. 

A. Flutter 


To provide a check for a special case of the pres- 
ent formulation and computer code, an aeroelastlc sta- 
bility analysis Is conducted for a tuned rotor with 12 
blades which was considered in Ref. (4). The blades 
have significant structural coupling Hue to the pres- 
ence of midspan shrouds. In addition to providing a 
check, the results give an example of the effect of 
structural coupling due to midspan shrouds on the sys- 
tem eigenvalues. The formulatiun in Ref. (4) Includes 
only single degree-of -freedom torsional blaHe motion. 

As a result, the bending frequencies in the present 
formulation are set to 10 times the torsional fre- 
quencies. This minimizes the inertial coupling be- 
tween the bending and torsional motions and provides a 
direct comparison with Ref. (4) results. The para- 
meters of interest for this rotor are given in Table I. 
The dependency of blade torsional frequencies on nodal 
diameter as taken from Ref. (4) is given in Table II 
for the reader's convenience. Figure 2 presents the 
eigenvalues for two cases of the tuned system; (1) no 
structural coupling and where the uncoupled blade tor- 
sional frequencies are equal to uq and (2) a flex- 
ible shroud ring in which the the frequency dependence 
is given in Table II. In this figure the plus symbol 
represents the eigenvalue for the case in which there 
is no structural coupling. The number adjacent to the 
symbol gives the predominant interblade angle for that 
eigenvalue. The arrow represents the location to which 
the eigenvalue moves (the end of the arrowhead) when 
the effects of structural coupling are considered. 

Note that for some modes the movement is so small that 
the arrows are not drawn. Within plotting accuracy, 
the results obtained using the present formulation are 
identical to those of Fig. 6 in Ref. (A). Note that 
there are different eigenvalue definitions between this 


paper and the reference. Since u only changes 
slightly, the addition of structural coupling has lit- 
tle effect on the stability of the Individual modes. 

As expected, for certain modes the frequency Is sig- 
nificantly affected. As was stated In Ref. (4), it is 
believed that this behavior Is representative of 
shrouded blades when the shroud ring is relatively 
stiff compared to the blades. 

To determine the effects of structural coupling on 
the coupling between the Interblade phase angle modes 
in the presence of mistunlng, this rotor was analyzed 
with 5 percent alternating mistunlng. With this type 
of mistunlng the torsional frequency alternates from 
blade to blade as you proceed around the disk. All 
odd numbered blades would have the identical frequency 
and all even numbered blades have the sar.te frequency. 
The percent mistunlng being defined as the difference 
in these two frequencies divided by the average. The 
resulting eigenvalues both with and without structural 
coupling are given in Fig. 3. The structural coupling 
generally only affects the frequency of the modes. 
However, it is interesting to note that the predomi- 
nantly 60 degree interblade phase angle mode (2 nodal 
diameter) is somewhat destabilized. 

To determine the effects of structural coupling on 
a more realistic rotor, an advanced fan stage called 
Rotor I is now analyzed. The parameters of this rotor 
are given in Table I. The assumed dependency of the 
tuned torsional frequencies (in vacuum) versus the no- 
dal diameter Is shown in Fig. 4. This assumption' is 
based on previous experience using finite element 
techniques to predict the natural frequencies of simi- 
lar bladed-discs. Figure 5 shows the effect of the 
structural coupling on the tuned system. As In the 
previous case, the effect is to significantly change 
the frequency of certain modes with little effect on 
the stability. This behavior is to be expected since 
each of the nodal diameter modes is uncoupled from the 
others. Note that there are still a number of un- 
stable modes. 

Kotor i is now analyzed with the assumption of al- 
ternate mistunlng. As seen in Fig. 6, all of the 
modes of Rotor I are stable in the presence of 10 per- 
cent alternating mistunlng without structural coupling. 
However, certain modes are significantly destabilized 
by the inclusion of structural coupling. Although all 
of the modes are stable, this behavior could result in 
a significant Increase In aeroelastlc response since 
the effective damping is decreased. 

Random mistuning Is now considered. Using the 
blade frequency distribution shown In Fig. 7 (mean of 
1.001 and standard deviation of 0.034, a typical pro- 
duction distribution) the eigenvectors of Rotor I with- 
out structural coupling were calculated and are dis- 
played in Fig. 8. The analysis was repeated including 
the effects of structural coupling and the correspond- 
ing results are shown in Fig. 9. Comparing Figs. 8 
and 9, it is seen that both the frequency and stabil- 
ity of the modes is significantly affected. It is 
important to note that the predominantly 45 degree in- 
terblade phase angle mode is moved into the unstable 
region. Although not shown, the blade amplitude dis- 
tributions for most of the modes are drastically dif- 
ferent when comparing the results with and without 
structural coupling.’ 

B. Aeroelastlc Response 

In the present formulation it is possible to con- 
sider an excitation from symmetrically spaced obstruc- 
tions located upstream from the blades. The number of 
obstructions is known as the engine order of the exci- 
tation and is usually represented by an integer f o 1 - 
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lowed by "E". For example, 5c would represent a per- 
fectly sinusoidal excitation resulting from 5 upstream 
obstructions. For a more detailed description of the 
formulation of the aeroelastic response problem see 
Refs. (1) and (2). 

To Illustrate the affect of structural coupling on 
response. Rotor I is again analyzed with the torsional 
frequency dependency shown in Fig. 4. The 5E resonance 
is investigated for the tuned, 10 percent alternating 
mlstuned, and randomly mistuned conditions. The Mach 
number is 0.8 and the reduced frequency is 0.6. All 
other parameters are the same as those given in 
Table I. These parameters were chosen so that all 
modes of the rotor are stable in the tuned condition. 

The eigenvalues for the tuned condition including 
structural coupling are shown in Fig. 10. The move- 
ment of the eigenvalue of Interest (-56 degrees inter- 
blade phase angle or 5 nodal diameter backward travel- 
ing wave) is represented by the arrow in Fig. 10. That 
is, the tail of the vector represents the position of 
the eigenvalue without structural coupling. The reso- 
nance peaks both with and without structural coupling 
are shown in Fig. 11. As expected, the frequency at 
which the rotor has maximum amplitude Is approximately 
4 percent lower with structural coupling. In addition, 
the maximum amplitudes are aporoximately equal. 

The eigenvalues for the system with 10 percent al- 
ternating mistuning are shown In Fig. 12. For this 
type of mistuning the 5E forcing function can excite 
only the -56 and 123 degree interblaot phase angle 
modes. The movement of these modes, depicted by the 
arrows in Fig. 12, represents primarily a change in 
frequency. As is the case for alternate mistuning, the 
single resonance peak is replaced by twin peaks which 
are shown in Fig. 13. Although not shown, the curves 
without structural coupling have the same general 
shapes and amplitudes with only shifts in frequency. 

The final case to oe investigated is that of random 
mistuning. The in vacuum frequency distribution of 
Fig. 7 is again used. The system eigenvalues are shown 
in Fig. 14. Comparing these eigenvalues with those 
without structural coupling (not shown), it is found 
that both the frequency and stability of certain modes 
was significantly affected. However, since the 5E 
forcing function excites all of the Interblade phase 
angle modes in this case, it is difficult to predict 
whether the amplitude of response will be larger or 
smaller than those of the rotor with no structural 
coupling. The response curve with multiple resonance 
peaks is shown Fig. 15. As expected, blade 2 has the 
maximum response at low forcing frequencies and blade 
17 has the maximum response at high forcing frequen- 
cies. Comparing this curve with that for no structural 
coupling (not shown), it was found that the inclusion 
of structural coupling had not significantly changed 
the overall response amplitudes. However, the blade- 
to-blade amplitude distributions have changed somewhat. 
For example, blades 10 and 14 have the highest ampli- 
tudes when the forcing frequency is approximately in 
the range 1 to 1.05. When there was no structural 
coupling, blades 3 and 13 were found to have the maxi- 
mum amplitude in this forcing frequency range. 

III. CONCLUSIONS 

The analyses of the effects of mistuning on coupled 
bending torsion flutter and response which were deve- 
loped in Refs. (_1) and (2) have been extended to in- 
clude the effects of structural coupling between the 
blades. The following conclusions are reached based 
on the parametric studies presented herein for bladed- 
Jiscs in which the motion is dominated by the blades. 


1. The addition of structural coupling can affect 
both th? frequency and aeroelastic stability. 

2. The stability is significantly affected only 
when the system is mlstuned. The affect can be either 
beneficial or detrimental depending on the type of 
mistuning, the type of structural coupling, and the 
cascade parameters. 

3. The resonant frequencies can be significantly 
affected by structural coupling in both tuned and mls- 
tuned systems, however, the peak response is signifi- 
cantly affected only for the latter. 
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TABLE I. - ROTOR PARAMETERS 


’ Parameter 

Symbol 

Ref. (4) rotor 

Rotor I 

Reduced frequency 

k 

0.75 

0.441 

Number of blades 

N 

12 

32 

Stagger angle 

c 

60* 

62* 

Elastic axis 
position 

a 

0.2 

0 

Gap-to-chord ratio 

s/c 

0.8 

0.777 

Mass ratio 

W S 

73.9 

141 1 

Radius of gyration 

r 

0.611 

0.47 

Elastic axis-C.G. 
offset 

X 

“s 

-0.2 

0 

Mach number 

M 

0.84 

1.18 

Bending to torsion 
frequency ratio 
tuned 

“hs'V 

10 

0.3 


TABLE II. - FREOUENCY 
DISTRIBUTION 
{IN VACUUM) 


r 

Tar 

0 

1 and 11 

2 and 10 

3 and 9 

4 and 8 

5 and 7 

6 

1.0000 

.99928 

.99724 

.98396 

1.0027 

1.0004 

1.0000 
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Figure 2. - Eigenvalues of rotor with 12 blades (tuned), 
M -0.84 
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Figure 3. - Eigenvalues of rotor with 12blades (5% alter- 
nating mistuning). M • 0.84. 



Figure 4. - Assumed frequency distribution of rotor I 
lin vacuum!. 
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Figure 5. - Eigenvalues of rotor I (tuned), M ■ 1. 18. 



Figure 6. - Eigenvalues of rotor I (10% alternating mis- 
tuning), M * 1.18. 
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Figure?. - Eigenvalues of rotor I. no structural coup- 
ling (random misfuning). M • 1.18. 
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Figure 11. - Tuned 5E response of rotor I. M ■ 0.8. 
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Figure 12. - Eigenvalues of rotor I. M ■ 0. a (10% alter- 
nating mistuning). 
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Figure 13. - 5E response of rotor I with 10% alternating 
mistuning, M * 0.8. 
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Figure 15. - 5E response of rotor I with random mistun- 
ing, M » 0.8. 


